The nature and boundary of the floating phase in a dissipative Josephson junction 
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We study the nature of correlations within, and the transition into, the floating phase of dissipative 
Josephson junction arrays. Order parameter correlations in this phase are long-ranged in time, but 
only short-ranged in space. A perturbative RG analysis shows that, in arbitrary spatial dimension, 
the transition is controlled by a continuous locus of critical fixed points determined entirely by the 
local topology of the lattice. This may be the most natural example of a line of critical points 
existing in arbitrary dimensions. 
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I. INTRODUCTION 

Systems where phases and phase transitions specific to 
a particular dimension are hidden in a higher dimensional 
manifold have recently come under renewed focusiiiS*2i4 
Of special importance is the case where the hidden phases 
and transitions are governed by spatially local, or zero- 
dimensional, modes4i^iSil A higher dimensional manifold 
with such local systems embedded in it may retain some 
of the properties that are effectively zero-dimensional. 

It has been shown recently that this is indeed the 
case in a dissipative Josephson junction array in any 
dimensioni-i^ For a system of Josephson-coupled super- 
conducting grains in a metallic substrate, the dissipa- 
tion is modelled by shunting resistors across the Joseph- 
son junctions.® The phase diagram of this model in the 
dissipation (Q;)-nearest-neighbor Josephson coupling [V) 
plane shows a "floating" phase at small a and V in which 
V is "irrelevant" in the renormalization group sense, so 
the phases of the superconducting order parameters on 
different grains are quantum disordered and the system 
is metallic.^ For larger a or V, V grows to strong cou- 
pling and the system orders into a superconductor. This 
phase diagram is shown in Fig. ^ 

In addition, it has been shown recently that all longer- 
range Josephson couplings, that may couple the spatially 
separated junctions in the array, are irrelevant in the 
metallic phase as long as V is irrelevant^ However, they 
all become relevant simultaneously with V at the same 
value of dissipation, etc, so that there are no further tran- 
sitions within the ordered phase; the ordered phase is a 
global superconductor. It has also been indicated that 
in the disordered metallic phase, the superconducting 
grains have long range in time, but only short range in 
space, correlations at zero temperature^ The metallic 
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FIG. 1: Phase diagram and perturbative RG flows for ID, 
2D triangular lattice (3D face-centered-cubic), and 2D square 
lattice (3D simple cubic). M and SC indicate metal and su- 
perconductor, respectively. The thick lines correspond to line 
of fixed points. Crosses indicate special points on the a = 
axis denoting {D + 1) -dimensional XY transition. 



phase has been called 'floating', the long-time-correlated 
zero-dimensional systems on the grains 'float' over each 
other. 

In the present paper we build on our earlier work4 by 
explicitly calculating the important correlation functions 
in the floating phase. By carefully including all relevant 
and dangerously irrelevant variables, we show that the 
long-time correlation functions of the order parameters 
on individual grains are indeed finite, whereas the spa- 
tially long-range correlation functions are exactly zero at 
zero temperature. At nonzero temperatures, the spatial 
correlation functions are exponentially decaying with a 
finite correlation length. Thus, each grain constitutes 
a dynamic local system on its own, and they are only 
short-range-correlated in space. It is then clear that the 
dissipation-driven quantum phase transition (QPT) be- 
tween the metallic floating phase and the global super- 
conductor is governed by the same (0 -I- l)-dimensional 
physics as in the 'quantum to classical' phase transi- 
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tion studied in the context of a single quantum parti- 
cle in a double wellSiiS or cosine-periodioii^ potential. 
Above a = Uc, the quantum fluctuations of the individual 
Josephson phases are effectively quenched by dissipation. 
Hence, in the absence of any quantum fluctuations, any 
infinitesimal V orders the extended system into a global 
superconductor. 

We further study this phase transition by computing 
the relevant renormalization group (RG) equations up to 
third order in V . While linear order recursion relations 
for V and a only show the existence of the critical ckc 
to separate the metal from the superconductor j^iSiii the 
higher order corrections reveal a continuous locus of crit- 
ical points in the {V — a) plane. This locus (see Fig. 
separates the metal {y irrelevant) from the superconduc- 
tor {y relevant), determining a phase boundary on the 
{y — ol) plane that approaches the F-axis for the spatial 
dimension D > \. The bending of the phase bound- 
ary toward the T^-axis is also what is physically expected 
and experimentally observedi^^ The functional form of 
the locus is determined entirely by the local topology of 
the lattice and does not depend on the global dimension. 
If the number of sides of the minimum closed loop on 
the lattice is three, for example in the 21? triangular or 
ZD face-centered-cubic lattice, the locus is a straight line 
making an angle with the a = line close to = 0. For 
the number of sides in the minimum closed loop greater 
than three, for example in ID square or ZD simple cubic 
lattice, the locus is a parabola close to y = 0. Since the 
forms of these loci do not depend on the global dimension 
in any way, they exist in any D\ this is a demonstration 
of a locus of critical points in any dimension; there is 
no upper-critical-dimension for this problem. These re- 
sults point to the intriguing possibility that the global 
transition in an extended system is ultimately driven by 
local physics at the level of a single junction. For Z? = 1, 
the line of critical points remains vertical, parallel to the 
a = line. This is analogous to the corresponding result 
in the zero-dimensional problem. Apart from intrinsic 
theoretical interest, these results also help unravel the 
nature of the floating phase and experimental observa- 
tions thereof. Local quantum critical points and their 
implications have also been widely discussed in the con- 
text of the cuprate high-temperature superconductorsi^ 
and heavy-fermion materials^ 

In section m we will give the preliminaries by intro- 
ducing the model and sketching the first order recursion 
relation for V . This will enable us to quickly identify the 
metallic floating phase as the one with a < ac, and the 
superconductor with oi> (Xc- In section ITTll we will dis- 
cuss the correlation functions in the floating phase. Some 
of the details are relegated to appendix IXI In section lTVI 
we will compute the corrections to the linear order recur- 
sion relation for V and derive the forms of the loci of criti- 
cal points for the triangular and square lattices in D = 2. 
Together with appendix^ which addresses the contribu- 
tion to the flows from third order cumulant expansion of 
the partition function, this completes the analysis of the 



perturbative renormalization group up to third order in 
V for the problem. The results are summarized in Fig.^ 
which also determines the phase boundary. In sectionlVl 
we will discuss possible experimental consequences of our 
results. Finally, the paper is summarized and concluded 
in section IVn 



II. THE MODEL 

We start with the following action describing a Joseph- 
son junction array coupled to dissipation in any dimen- 
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(1) 



Here the sum is over nearest-neighbor pairs, and 

we will refer to A9ij ~ {0i — 0j) as the phase difference 
on the bond A0y (w„) is a Fourier component of 

A^y (r), and a;„ = ± 27rn//3, with n an integer and f3 the 
inverse-temperature. C, the capacitance of the supercon- 
ducting grains, is assumed to be small, and so Eq ^ 1/C 
constitutes the largest energy scale in the problem. The 
dimensionless variable a = Rq/R, where R is the shunt 
resistance and Rq = h/ie^ is the quantum of resistance, 
couples the Josephson-phases A9ij to a local dissipative 
heat bath, with a Ohmic dissipation.* 

To establish the existence of a dissipation-tuned quan- 
tum phase transition (QPT), we employ a frequency-shell 
RG which is perturbative in V. Let's divide the held Oi{T) 
into slow and fast components 9is{T) and ^^/(t), such 
that e^r) = ;^E|.„|<../b^~^K)e'""^ and %(t) is 
given by a similar expression where the sum is between 
LUc/b < < u>c- Here, uJc ^ Eq is the high-frequency 
cut-off of the frequency integrals, and 6 > 1 is a frequency 
rescaling factor. 

We can now write the partition function Z as a func- 
tional integral. 



Z = N 



-S/h 



-Sq /h-\-\n{e 



(2) 



Here N and N' are normalization constants, Sq is the 
slow-frequency component of the quadratic part of the 
action containing the flrst and the third terms in Eq. 
S' contains the second term, and {■■■)of denotes average 
with respect to the fast component of the quadratic part. 
After computing the averages, we rescale r, t' = r/b, 
to restore the original frequency cut-off, and then rede- 
flne the coupling-constants to complete the renormaliza- 
tion. The dissipation term is dimensionless, and so it 
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is held fixed in RG. The first term in Eq. (QJ then has 
T-dimension —1. Since it is not renormahzed by any of 
the other terms, the couphng-constant C renormahzes to 
zero by power counting. As we will show in the next sec- 
tion, we can put this term equal to zero for the purposes 
of the RG, since inclusion of it does not change any of the 
results. For computations of some correlation functions 
in the floating phase, however, C is dangerously irrele- 
vant; some of these functions will depend singularly on 
C. 

The average in Eq. ^ can be easily performed in the 
leading order in V, 

(exp[y^ / rfTcosA6iy(T)])o/ 

= exp[V fdr Y: e-^^iM(e-^«i(r))o/] 

= exp[y V / dr COS A0i J {T)e 
,: Jo 



-i(Ae{,.(r)^)o/l 



where, 



(AC(r)^)o/ 
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(4) 



Here zq — z/2, where z is the coordination number of 
the lattice. In deriving Eq.^jwe have taken C to be zero, 
which, as mentioned before, is allowed for the present 
purpose. After rescaling the time, r' = t/&, and taking 
6 = e', where ^ > is infinitesimal, we end up with a 
linear order RG flow equation for V , 



al zoa 



(5) 



It is clear from this equation that a — Uc — I/zq is 
a critical fixed point for the flow of V; for a < ac, V 
scales to zero, while for a > V grows. In the phase 
where V scales to zero, the junction-phases are quantum- 
disordered and the system is metallic by construction due 
to the existence of the shunt resistors. When V flows to 
higher values, the system is phase-ordered and supercon- 
ducting. 

Now we show that in the metallic phase, all longer- 
ranged Josephson couplings beyond the nearest-neighbor 
coupling V, which can couple spatially separated A^^j's 
on the array if included in the starting action, are irrele- 
vant. These couplings are described by the general term. 



1 Jo 



Here, 9{fi, r)'s are the order parameter phases in the su- 
perconducting grains, the Si's are arbitrary vectors sep- 
arating the lattice points on the lattice and Si is an 
integer-valued function of the layer number i satisfying 
Si = 0. The last constraint ensures the absence of an 
external fleld, which implies "rotation invariance" under 
adding the same constant to all of the 9's. Note that the 
special case sq = +1, Sq = 0, si = —1, Si = a^, where 
is a nearest neighbor vector, is just the nearest neighbor 
Josephson couphng V in Eq. 

Calculations analogous to that leading to Eq. jS)) give, 
to first order in J, 



dJ{[s,]) 
dl 



r(N). 



with 



(7) 



(8) 



The "potential" U{f) = ^ - 1)//,-, with = 

^^(1 — e**"T), where the sum over 7 is over all nearest 
neighbors. It is straightforward to show that U{r) is 
the "lattice Coulomb potential" of a unit negative charge 
at the origin, with the zero of the potential set at r = 
0. That is, U{r) satisfies the "lattice Poisson equation": 



zU{r) = — <5^g. For a symmetrical (e.g.. 



square, hexagonal, cubic) lattice, the left-hand side is just 
the "lattice Laplacian", approaching V^C/ x 0{a?) where 
(a EE \d^\). 

The quantity F in (TJ and © is then just equal to the 
potential energy of a neutral (since Si = 0) plasma 
of quantized (since the siS are integers) charges Si on 
the lattice. The most relevant J([si]) is clearly the one 
that corresponds, in this Coulomb analogy, to the lowest 
interaction energy. Note that strictly speaking F cor- 
responds to twice this energy, because the sum in (jSJ 
double counts. Apart from the trivial configuration in 
which all the Si = , the lowest energy configuration 
is clearly the one in which there are two equal and op- 
posite unit magnitude charges on nearest neighbor sites: 



i.e.. So = +1, 1^0 = 0, si 



As we discussed 



(6) 



earlier, this corresponds to the nearest-neighbor Joseph- 
son coupling in equation Thus, it is established that 
that coupling is, indeed, the most relevant, as we asserted 
earlier. Furthermore, using simple symmetry arguments, 
we can show that for a symmetric lattice (e.g., square, 
hexagonal, cubic), where all nearest-neighbor sites are 
equivalent, U{dj) — —-^i which recovers the recursion 
relation (0) for V . As a result, all other couplings are 
irrelevant for a < = '^/zq; hence, they affect neither 
the floating phase nor the transition between it and the 
{D + l)-dimensionally coupled phase. 

Thus, the metallic phase is in a sense spatially decou- 
pled. This decoupling of space, but finite correlation in 
time, can be more clearly seen by the correlation func- 
tions in the floating phase discussed in the next section. 
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III. CORRELATION FUNCTIONS IN THE 
FLOATING PHASE 

In this section, we first show that the variable C is dan- 
gerously irrelevant in Z? < 2 for some correlation func- 
tions in the floating phase, but not for the ones required 
for the RG. Including all relevant and dangerously irrel- 
evant variables, we will then compute various long-time 
correlation functions fEas. 1131 1141 1151 and show that 
they are flnite. Finally, the spatially long-range corre- 
lation functions (Eqs. ^1 1211 I23|) will be shown to be 
exponentially decaying with finite correlation lengths at 
a non-zero temperature and exactly zero in the limit of 
zero temperature. Thus, the intrinsically local and dy- 
namic character of the floating phase will be established, 
indicating that the dissipation-driven QPT between this 
phase and the global superconductor is governed by the 
same (0 + l)-dimensional physics as in the corresponding 
zero-dimensional transitionn^iifliiiiiS, 

From the gaussian part of the action in Eq. ^ we find 
that 



1 



Cc.2 + (a/27r)/(k)K 



(9) 



Here /(k) = |(e*'' * — 1)^, where 5's constitute the set 
of nearest neighbor vectors. Since /(k) ^ zok^ for small 
k, it is easy to see that in dimension D < 2, wavenum- 
ber integration of the right-hand-side (RHS) of Eq. |51di- 
verges in the infrared for C ^ 0. This divergence forces 
us to include a non-zero C for the correct evaluation of 
some of the correlation functions in the floating phase, 
even though C is an irrelevant variable in the RG sense. 
The recursion relations, however, remain unaffected as 
we show below. 

The only correlation functions we need for the pur- 
poses of the RG involve the variables defined on a bond, 
A9ij{T), the simplest of which is given in Eq. 01 where 
we have taken C to be zero. In calculations higher 
order in V, we will need higher order correlators 
{A9ij{T)A9ki{T))Qf, still among the bond- variables. For 
such functions, it is easy to see that including a non-zero 
C on the RHS of Eq. |51 only amounts to introducing a 
high-frequency cut-off Eq ^ 1/C in the frequency inte- 
grals. In the case of {Adij{T)'^), summing over j over the 
nearest-neighbor vectors assuming symmetry among the 
nearest neighbors, dividing by zq, and using Eq. IHl we 
get in D = 1, 



1 rik 

/(k) 



1 f'^ dk 



E 



zoJo 2^ 1^ Cc.2 + (a/27r)/(k)K|' 

where A is a momentum cut-off. It is easy to see that 
since /(k) appears both in the numerator and the de- 
nominator, the infrared divergence in the wavenumber 



integration is eliminated. Since in the limit of small k the 
integral is zero, /(k) is order one in the integral, hence 
Cw^ can be neglected with respect to (Q;/27r)/(k)|w„| 
for small uon. An explicit evaluation of the integral 
shows the same logarithmically divergent behavior as in 
Eq.l (A0y (T)2) = -l-ln(/3So), where Eq - K/C and 

\nK = 1^ ln(Q;/(k)). Hence, for the purposes of the 
RG, taking C = in the expression for < |6'(k, LJn)^ > 
and taking a compensatory high-frequency cut-off ~ 1/C 
are indeed justified. 

From {AOijirY) ^ we can calculate 



(exp(igA%(r))) -l/(/3So)^ 



(10) 



Here, g'^ ="gz (? ~ '^)^' ^^"^ g is a real number. Note 
that for g = 1, the simplest such function that one can 
construct, the RHS of Eq. ^| does not vanish at zero 
temperature since the corresponding q' vanishes. This 
is analogous to a single paramagnetic spin developing a 
non-zero expectation value in the presence of an applied 
external field. In the present case, the Josephson cou- 
pling V acts as the field inducing a non-zero expectation 
value of e*^^*^ , even in the complete absence of all other 
couplings. This expectation value does not constitute an 
order parameter since the broken symmetry is not spon- 
taneous, but merely induced by the applied field. It forces 
us, however, to introduce the parameter q to get corre- 
lation functions which do decay algebraicalljii as T ^ 0. 
Note that q'^ is non-zero as long as q is not an integer. 

The result for q' is analogous to a well-known result 
in surface rougheningfiii^ which we shall sketch here 
for the sake of completeness. The result, (A0y (r)^) = 
2\n{[3EQ)/ z^a, correctly gives the low-temperature fluc- 
tuations of AOij throughout the floating phase, up to an 
unimportant additive constant. Since this is valid irre- 
spective of the actual bare value of V ^ we can apply it 
even in the limit of large V . But since in this limit AOij 
is quantized in integral multiples of 2tt, AOij = 27m with 
n an integer, we can write for the function in Eq. 1101 

oo 

(exp[igA6'ij(T)]) = ^ P{n)exp{27Tiqn) 



n— — oo 



j: ( 



TTZQa 

Vln(/3So; 



exp 



2'Kiqn 



HPEq) 



(11) 



Here P{n) is the gaussian distribution function for n with 
(n) = ^(A0ij(r)) — 0, the standard deviation (7„ — 
^(A%(r)2), and the sum over n, instead of an integral 
over a continuous variable, takes care of the quantization 
condition on AOij. Doing this sum using the Poisson 
summation formula gives 

oo 

(exp[zgA%(r)]) ^ (/3i?o)-('-^)'/(^""). (12) 

S — — OQ.S^Z 

Clearly, the dominant term in this sum as T ^ is the 
one with the smallest value of (g — s)^. Keeping just this 
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term yields Eq. H10|l . This shows that this correlation 
function indeed goes to zero as the temperature goes to 
zero. The unequal-time correlation function of the bond- 
variable is also algebraic, 



(exp[z(7(A%(T)-A%(0))]) 



1 



(£;oT)29'V^oa ■ 



(13) 



The correlation function of the bond variables alge- 
braically decaying in time, Eq. 1131 is valid in all dimen- 
sions. The same is not true, however, with the correlation 
functions of the site variables, ^^(t), where i is a grain 
index, since now for D < 2 the wavenumber integration 
of < |6'(k, w„)p > is divergent as C — > 0. For D > 2, 
the effect coming from C is subleading and the long-time 
correlations remain algebraically decaying. In D = 3, for 
the spatially local but unequal in time correlator we get, 
ignoring unimportant constant factors. 



(exp[ig(0,(r)-0,(O))]) 



1 



,'2a 



(14) 



where A is the momentum cut-off. In D 
correlation function is given by. 



2, the same 



{exp[iq{e^{T) - e^{0))]) ^ exp 
while in D = 1 it is, 

(exp[zg(0,(T)-0,(O))]) -exp 



(in * 







2nC 



(15) 



(16) 

Even though the long-time behavior is algebraic only 
for the bond- variables, Eq. El and for the site vari- 
ables only for D = 3 and above, Eq. ^1 temporally far- 
separated sites are still correlated for Z3 < 2 according 
to Eqs. ifTsl [TBI) . However, long-range spatial correla- 
tion functions of the grains, (exp[i(7(6'j(T) — Sj+r ('''))]) j 
are identically zero at zero temperature in any dimen- 
sion. Spatially far-separated grains are completely un- 
correlated, each one of them constitutes a dynamic local 
system on its own. In this sense, the temporal and spa- 
tial behaviors of the system are completely decoupled at 
T = 0. It is also a manifestation of the intrinsic local na- 
ture of the problem. At nonzero T, the spatial correlation 



functions are exponentially decaying in all dimensions; 
the grains are only short-range-correlated. 

As we show in appendix ^ the correlation function of 
the phases of the order parameters on the grains follow 
the simple scaling law. 



(0,(r)0,+,,(T)> 



27rr 



2-D 



azo 



(17) 



where, ^ is the correlation length at nonzero T, ^ — 
(azo/27rCT)^/2. As we will see below, even though f 
diverges as T 0, there is very little correlation left 
between the order parameters themselves at low temper- 
atures, and eventually the order parameter correlation 
function vanishes at T = 0. In Z? = 3, we find. 



9ir/0 



27r2 r 
1 

2^ 



2 „ exp(-r/^). 



r << ^ 



(18) 



This gives for the correlation function of the order pa- 
rameters, 



— -e « , r » ^ 

nazQr r 



— - 1 



, r « ^ 
(19) 



where the subscript c denotes the connected piece of the 
correlation function. It is clear that as T ^ 0, these 
functions vanish for large r. Hence, there is no corre- 
lation among the spatially far separated grains at zero 
temperature. In _D = 2, we get for the scaling function 



1 f 3 

-— ^(-)3exp(-r/0, r»e 
V(27r'^) r 



27r2 



[ln(2^/r)]^ 



(20) 



The order parameter correlation function becomes. 



(exp[i9(0j(T) - ei+r{T))])c 



exp 



exp 



A'kuzq 



TT azo r 4 







(exp[ 



{q'Hr /20f 
nazQ 



r << ^ 



(21) 



r 



It is straightforward to check that this function is also zero at T = 0. Finally, in D — I, g{r/^) is given by. 



gir/0 - (l/7r)(e/r)2exp(-r/C), 



r >> ^ 



--ln(-) + -^exp(-r/0, r « ^ (22) 
TT r TT r 



6 



The corresponding order parameter correlation function is, 



(exp[z(7(0,(r)-0,+,(r))]). 



exp 



^/^)(2g'V)/(a2;o) _ 

I 



(^)(L)exp(-r/0, 
azQ r 



exp 



r » £, 
r « ^ 



(23) 



This function also approaches zero as T ^ 0. Thus, we 
have established that the spatially far-separated grains 
are uncorrelated at T = 0, and correlated only with ex- 
ponentially decaying correlations at a non-zero tempera- 
ture. 



IV. RG AT SECOND ORDER IN V 

To establish the existence of a locus of critical points 
in the a~V plane, we need to compute the higher-order 
corrections, if any, to the first order recursion relations 
for V, Eq. (O. Performing the cumulant expansion of 
Eq. (|2Jl in second order in V, we get 



exp 



2 



cos A0ki{T'))af 



dr^ COS Ae,j{T))lf 
+ A0^;(r'))(e-<^<(")^''^'("')>«^ - 1) + cos(A0?.(r) 



exp 



(24) 



To avoid the generation of spurious long-ranged be- 
havior in the correlation function of the fast modes, 
{A9{j{T)A9li{T'))of, we need to adopt a smooth cut- 
off prescription^ in our frequency-shell RG. For the 
unequal-time correlation function on the same bond, this 
can be written as 

Gf,., ((r-r'),5)^(A0/,(r)A0/,(r'))o/ 



— y 

|i^„l<ti. 



exp(zu;„ T - r //^c^w„^^ 
jTTi \ / '^"/ -r))^25) 

a/2TT\uJn\ b 



where f{x) is a smoothing function with the properties 
f{x) ^ for a; << 1, and f{x) ~ 1 for a; >> 1. The 
choice of f{x) is somewhat arbitrary as long as it entails 
the properties, Cf^- j^- ((r— r'), b) ^ for |r— t'| oo, and 

G{j y ((t — r'), 6) ^ for & ^ 1, to the correlation func- 
tion of the fast modes. A specific choice is given in Ref. 



Ha . where G goes to zero exponentially with |r — r'| for 
large |r-r'|, and G{,- ^ (0, b) = This does not mod- 

ify any of our previous calculations and we will assume 
that G{j fci((''' — t'), 6) is short-ranged in |t — r'|. The t 
and t' in the integrals of Eq. H24|) are now constrained to 
be close to each other and we can use gradient-expansion 
in (r — r')to simplify the terms within parenthesis. It is 
clear that the second term does not renormalize V. The 
first term, however, upon gradient expansion generates 
cos(A6'|j(t) + A91i{t)) summed over all bonds {i,j) and 
{k,l). 

It is clear that there is no renormalization of V in 
second order in one dimension. In the sum over 
and {k,l), if (j,j) — {k,l), one generates a term ~ 
cos(2A6'fj(r)) which is irrelevant close to Uc- The second 

term also produces an irrelevant term, ~ {-—g^r-^Y ■ For 
(*) j) 7^ Oj too, it is easy to see that there is no renor- 
malization of V, since cos(A6'y + A9ki) does not produce 
cos(Af?m„); in one dimension there is no closed loop on 
the lattice to do this. Contrast this with D = 2 triangular 
lattice, where because of the existence of closed loops in 
the lattice cos(A0m„) is produced from cos(A6'y + A9ki) 
when the bonds {i, j), {k,l) , {m,n) form the three sides 
of a minimum triangle. Similar effects exist for any lat- 
tices in D > 1. Coming back to D = 1, we have checked 
explicitly, as shown in the appendix ^ in a more general 
setting, that there is no renormalization of V also at the 
third order. Thus we speculate, and confirm up to third 
order expansion in V, that in _D = 1 there is no correction 
at all at any higher order to the linear order recursion re- 
lation Eq. El This was conjectured in the context of the 
corresponding 0-dimensional problem in Ref^, here we 
have extended it to one dimension. Notice that in Eq. ^ 
because of the existence of the C-term, the one dimen- 
sional problem does not map onto the zero dimensional 
problem in any obvious way. We will see below that for 
D > 1 there are corrections to Eq. ||SJ) at higher orders. 

In 2D triangular array, and in fact in arbitrary dimen- 
sions where the minimum closed loop is a triangle, there 
is a contribution at order V'^ to the recursion relation 
for V. This comes from the first term in (|24ll . where, 
in the sum over {k,l), when (fc,/) is the bond adjacent 
to (i,j) in a minimum triangle, the sum over and 
(fc, I) produces a single sum over all bonds in the lattice. 
This is illustrated by Fig.(la). If =< 1,2 >, and 
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{k,l) =< 2,3 >, then cos(A6'fj- + A6l^,;) = cos(A6'^3) since 
+ A6'23 + A03i — around a triangle. The first term 
within parenthesis in Eq. 124() can then be written as 

e-'^Yl drY,cos{Aet,{T)) 
X / d(r'-r)[cxp(-G(2,23(r-r'))-l], (26) 



where 23 ^ t"') is the nearest-neighbor-bond cor- 
relation function of the fast modes. For the triangle 
in Fig. (la), using A6'i2 -I- A6'23 + A6'3i = 0, we get 

G(2,23(^ - r') = -(1/2)G{2,i2(^ - r'). The precise val- 
ues of these correlation functions are somewhat arbi- 
trary since they depend on the choice of the smoothing 
function f{x). Assuming they are finite only close to 
|r — r'l = 0, we can then define tq so that the (t' — t) 
integral in Eq. is given by To[exp(iG(2 12(0)) — 1]. 

and rescaling r, we find the re- 



Using 6(2,12(0) 



21nb 



cursion relation for V at second order in V 
al zoa 



(27) 



where Ci is a positive constant Ci = ro/4zoa. In 2D 




4 r 



Cb) 



FIG. 2: The minimum closed loop generating higher order 
renormalization of V for (a) 2D triangular lattice, and (b) 
2D square lattice. The next-nearest-neighbor interaction in a 
square lattice is denoted by J. 

square lattice, and also in a lattice in any dimension 
where the minimum closed loop has a number of sides 
greater than three, the first non-zero contribution above 
the linear order in V occurs at order V^. Interestingly, 
this can already be seen at the level of second order cu- 
mulant expansion of the interactions. This is illustrated 
by Fig. (lb). Note that in the square lattice, the first 
term in Eq. (|24|l . upon the usual gradient expansion, gen- 
erates a next-nearest-neighbor interaction J. In a consis- 
tent RG-treatnient, J (and in fact all such longer-ranged 
Josephson interactions) should be included in the start- 
ing action. As has been shown before, they are all ir- 
relevant in the floating phase as long as V is irrelevant. 
However, as soon as V becomes relevant at ac, their re- 
cursion relations pick up source terms of order . They 
in turn generate higher order corrections to the recursion 
relations for V. 



It is clear that the lowest order correction comes from 
the next-nearest-neighbor interaction J. Including this 
term in the starting action in Eq. J^l, and with the help 
of Fig. (lb) and Eq. (|7J), we find for the recursion relation 
for J, 



dl ^ 



1.26 

Zoa 



)J+i 



0.19r, 



ZQa 



(28) 



where Tq is defined via J^^d{T — T)[exp(— 23(''" ~ 

r')) — 1] = ro[exp(— 0(2 23(0)) ~M- By an exactly analo- 
gous calculation, we get for the recursion relation for V, 



at zqu ZQa 

where Tq is defined by a third integral, /_f d,{T — 
r)[exp(-G{3,34(T ~t'))-1] = To [exp(-G(3,34(0)) - 1]. 

dj 



Solving for ^ = at a = etc, we find J ^ V , and 



putting it back in Eq. (|29|l we get. 



dV_ 
dl 



(1 )V + C2V' 

Zoa 



(30) 



where G2 is a positive constant given in terms of the time 

cut-offs, G2 = ( °-^r"""' ). 

There is a possibility of an order term appearing in 
the recursion relation for V at third order expansion of 
the interactions for the square lattice. This contribution 
may come from a term of the form 



24 



rP rP f0 

/ dn / dT2 / dT,Y^J2Y. cos[A0^,(ri 

Jo Jo JQ I- A ,, n / \ 



+ A0^.,(r2) + A0^„(r3)]e-^ [exp(G4,,(Ti - r^) 

+ ,^„(t2 - T3) + Gf„„ ,^-(t3 - Tl)) - 1], 

which appears in the third order cumulant expansion of 
Eq. lI2Jl. We have checked explicitly, as shown in the ap- 
pendix^ that this docs not happen. In fact, the contri- 
bution to Eq. (|30|) from third order cumulant expansion 
is precisely zero. Thus there is no cancellation of the 
order term we have found above. Note that for any 
lattice, where the minimum closed loop has number of 
sides greater than three, the order term appears via 
the generation of the next-nearest-neighbor interaction 
in the RG. Thus 2D honeycomb lattice, for example, will 
also have a recursion relation for V similar to Eq. |I3C 



Combined with the flow equation for a, ^ = 0, which 
is true in our perturbative treatment close to a^ for all 
lattice types, we have sketched in Fig. the perturbative 
RG flows for ID, 2D triangular lattice and 2D square 
lattice. 
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V. EXPERIMENTAL CONSEQUENCES 

Certain aspects of our theory can be seen in a recent 
experimenlii^ (see also references therein) , although some 
caution must be exercised. In this paper, a phase tran- 
sition is shown to exist at Uc — ^/zq, for small values 
of V/Eq, where zq = 2 for the square lattice fabricated 
in the experiment. This is consistent with the predic- 
tions first given in Ref. and later in Ref. ^3 from 
a single junction approach. Moreover, the experimental 
phase diagram, especially the absence of a dissipation- 
driven transition above a certain value of V/Eq, bears 
a remarkable resemblance with what can be trivially in- 
ferred from our Fig. There are two issues that deserve 
further attention: (1) the quantum disordered state of 
the resistively shunted Josephson junction array is metal- 
lic, not an insulator, as inferred in this experiment. This 
is simply because the metallic shunts will short the cur- 
rent when the superconducting phase difference between 
the junctions become incoherent, unless disorder effects 
in two dimensions turn this metallic behavior to an in- 
sulating behavior. However, the experimental system is 
claimed to be free of disorder. This metallic behavior was 
not tracked down in this experiment, but only the insulat- 
ing behavior of the unshunted array; further experiments 
can clearly be revealing; (2) the correct model of the ca- 
pacitances must also include longer ranged contributions, 
as they are obviously present in the experiment. 

Given the novelty of the floating phase, further ex- 
periments will be a worthwhile effort. For values of 
R > Rq = , we expect the system to be in the floating 
phase, where the grains are only short-range-correlated 
in space. We expect the current-voltage characteristic in 
this case to be a non-universal power law controlled by 
dissipation a, similar to the results obtained previously^*^ 
in the context of isolated junctions. For R smaller than 
the quantum of resistance, the system will be in the 
superconducting state, and the temperature-dependent 
non-universal power law will be similar to the classical 
XY model due to vortex unbinding in the presence of an 
imposed currentuSi Recall that V, along with all other 
longer range Josephson couplings,^ flow to larger values 
in this regime. An experimentally testable prediction for 
a non-universal exponent v controlling the critical scal- 
ing of a characteristic temperature Tch, Tch oc {SV)" with 
6V = \V—Vc\ and Vc the critical value of V on the (V—a) 
plane, was made in our earlier work.** If verified, this will 
be a clear cut signature of the phase transition studied 
here. By varying the strengths of the Josephson coupling 
and the shunt resistors, the critical locus and hence the 
phase boundary on the (V — a) plane can also be directly 
examined. 

The ideas presented here assume significance also in 
a broader context In a bid to understand some puz- 
zling experimental observations in the context of heavy 
fermion materials, and the cuprate superconductors 
QPT's are sought where extended systems may go 
through global phase transitions which are ultimately 



driven by inherent local ones. The transition we have 
studied in the dissipative Josephson junction array pro- 
vides just such an example, albeit only in spirit, since we 
have not treated the fermions here. It will be interesting 
to extend the ideas and results found here to fermionic 
problems to make contact with a wide range of systems 
of experimental interest. 



VI. SUMMARY AND CONCLUSION 

In this paper, we have clarified in detail the nature 
of the correlation functions within, and the transition 
into, the floating phase of a dissipative Josephson junc- 
tion array. We have shown that the long-time correla- 
tions in the floating phase are flnite locally for each su- 
perconducting grain, but the spatially long-range corre- 
lations are identically zero at T = 0. At any nonzero 
temperature, the spatial correlation functions are expo- 
nentially decaying with a flnite correlation length, indi- 
cating locally dynamic grains which are only short-range 
correlated in space. These correlations help identify the 
dissipation-driven metal to superconductor transition in 
the extended system of the array with the 'quantum to 
classical' dynamic transition existing at the level of a sin- 
gle junctioujSiiS the local transition drives the global one. 

Further, we have shown that the transition in the ex- 
tended system is controlled by a locus of critical points in 
the a — V plane with continuously varying critical expo- 
nents. To derive our results, we computed the perturba- 
tive recursion relations for the nearest neighbor coupling 
V up to third order in V. Since V is the assumed small 
parameter, the calculation is valid only close to the crit- 
ical value of a ~ ac needed for ordering with inflnites- 
imal Josephson coupling. We have found no correction 
to the linear-order recursion relation, Eq. [S] for D — 1. 
In this case, then, the fixed line close to Uc is parallel 
to the V-axis. This is a result similar to that obtained 
for the corresponding zero-dimensional problem treated 
m Refd Here we have extended their result to D = 1. 

For any £) > 1, however, there are higher order correc- 
tions to Eq. Coupled with the recursion relation for 
a, ^ — 0, the locus of critical points bends toward the 
V^-axis determining a phase boundary. This locus is a 
straight line in any D > 1 as long as the number of sides 
in a minimum closed loop is three, it's a parabola if that 
number is greater than three. The functional form of the 
locus, which determines the phase boundary between the 
metal and the superconductor, is entirely determined by 
the local topology of the lattice, and not on the global di- 
mension. Consequently, the locus of critical points exists 
in any dimension; there is no upper-critical-dimension for 
the problem. All of these results are consistent with the 
global transition being driven by the inherent local one. 

The correlations in the fioating phase and the phase 
boundary studied here can be experimentally probed by 
measuring the current-voltage characteristics with the re- 
sistance of the shunt resistors continuously tuned. In a 
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broader context, taken together with our earher result^ 
that all longer-range Josephson interactions induced by 
the lattice among a set of periodically spaced local sys- 
tems, which each goes through individual local transi- 
tions with dissipation as the tuning parameter, are irrele- 
vant in the RG sense, our results here assume significance 
in the context of the search for local quantum criticality 
in extended systems i 
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APPENDIX A: SPATIAL CORRELATION 
FUNCTIONS IN THE FLOATING PHASE 

To compute the spatial correlation functions of the 
phase variables at a temperature ^, we need to evalu- 
ate the following integral in frequency and wavenumber 
space: 



d'^k exp(ik.r) 
{2nY a7(k) + C|c^„ 



(Al) 

where a' — ajl-K for notational convenience. Noticing 
that only k ^ \/r modes will contribute significantly to 
the integral, and for large values of r, we expand the func- 
tion /(k) ~ z^k'^ in the denominator. The fc-space inte- 
gral is then simply the Fourier-transform of the Yukawa 
potential at wavenumber space. 



1. Correlation function in D=3 

In _D = 3, doing the fc-space integral, which simply 
gives the three dimensional Yukawa potential, and after 
a change of variable, we get 



.{r)e^+r{T)) = 



1 1 



du 



exp(-u), (A2) 



27r2a'zo r J^^^ u 

where = C(T)^^ — The w-integral produces 

the exponential-integral functionj24 and the exact corre- 
lation function is given by. 



{e,{T)e.+r{T)) 



1 1 



27r2 



a' zq r 



m-r/0- 



(A3) 



Taking the asymptotic values of Ei{x) for large and small 
|a;|, we end up with Eq. ^| With the help of the result 



2. Correlation function in D=2 

In I? = 2, the angular integration in fc-space produces 
Jo(fcr), where Jo{z) is the Bessel function of the first kind 
of order zerop24 The remaining part of the fc-space inte- 
gral produces the function 2'KKo{{-^j^)^r), where Ko{z) 
is simply related to the Hankel function of the first kind 
Hq{z). After a change of variable, we end up with the 
integral, 

{e.{r)e,+rir)) = / -Ko{u). (A4) 

For r >> ^, we use the asymptotic expansion of the 
function Kq{u) for large u, and after doing the it-integral 
we get the result 



1 



V(;)(^)3/4exp(-r/20 



X W^-3/4,-l/4(''/0> 



(A5) 



where W\^f^{z) is the Whittaker function. Asymptoti- 
cally expanding the Whittaker function for r >> ^, we 
end up with the first line of Eq. [501 In the limit r << ^, 
we notice that the integral in Eq. IA4I is dominated by 
small u. Expanding Ko{u) for small u, and doing the 
M- integral we get the remaining part of Eq. [201 With the 
additional result in D 
we finally get Eq. ED 



2, {(^dr?)-s^ilni^?, 



3. Correlation function in D=l 

In D ~ 1, after doing the fc-space integral and a change 
of variable, we get, 



na 
r 



— \ — exp(-u) 

'Zo Jrii 



na zq 4 



(A6) 



Asymptotically expanding the exponential-integral func- 
tion for large r/^, the first term in the expansion cancels 
the second term within the square bracket of Eq. IA6I 
The first subleading term in the expansion produces the 
first line in Eq. Expansion of Ei{r/£) for small 

values of the argument produces the remaining part. 
With the additional result that in D ~ 1, < 0i{TY >— 

Jl/2 



(3 ' , we end up with Eq. 



Finally, by doing the change of variables (3ljj — Vl and 
kr = Q where r is the magnitude of the spatial coordi- 
nate r, and taking 1 / {a' zq) outside the integral in Eo. lAll 
we get the scaling relation Eq. El 
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APPENDIX B: RG FROM THE THIRD ORDER 
CUMULANT EXPANSION IN V 

Here we show that the term of order in the recursion 
relation for V , found from the cumulant expansion of 
the partition function of Eq. |21at third order, identicahy 
vanishes for any lattices in arbitrary dimensions. Using 



the formula for third order cumulant expansion, 



(exp[/])=exp[(/> + -((/^>-(/)^)(l-(/)) + ^-^, 

(Bl) 

We get, for / = y /o''rf^E(.,,) cosA%(t), 



exp[/] = exp l/e-'"^/^°" / dr ^ cos A0^(t) + (^e-^^^^/^oa ^ ^ / / cos(A0^,(ri) + A0^;(t2)) 



^0 



1) + cos(A0|^.(ri) - A0ii(T2)) X (eGf,..(n-r.) _ _ Ve-'-'/^'"^ J^t^ ^ cosA0^„(t3) 



(m,n) 



(/^>o/ {f)h 



(B2) 



After evaluating the correlation functions of the fast 
modes, G{jj,;(tiT2)'s, doing gradient expansions in {t2 — 
Ti) and (t3 — Ti) (it can be seen by straightforwardly 

writing out ( " ^ ^) in the above expression that, 

all terms where ti , T2 and T3 are not constrained to be 
close to each other cancel among themselves) , and finally 
rescaling ti , a general term that may contribute to renor- 
malize V at third order has the form (with ln6 = 

y3^l^~3l/zo^ r dT^Y^Y. E COs(A%(ti) 

+ AOkiiri) + Ae„,n{Ti)] X f d{T2 - ri) 

\ d{T3 - Ti)(exp[-G{^f^i{T2 - ri) - G{;,„„(t3 - T2) 
- GL,,(^3-ri)]-l) (B3) 

Note that the Gf^H (■^i'^2)'s are evaluated within a thin 
shell o^cG ' < 1^1 < around the high frequency cut- 



off. Using the property that they go to zero as Z ^ 0, we 
have to expand their exponentials appearing in the above 
expansion in linear order in the infinitesimal parameter I 
to get their contributions to the recursion relation of V. 
Hence, we do this expansion in / right away in all of the 
terms appearing in Eq. IB2I We also use the trigonomet- 
ric identity cos(A) cos(B) = i[cos(A + B) + cos(A - B)] 
to simplify the second term of the argument of the expo- 
nential in Eq.mi Then writing out (^^^^ - -^4^) using 
three variables €1,62, £3, which take the values ±1, per- 
forming the averages with respect to the fast modes (note 
that no higher order correlation functions arise, only the 
familiar Gf^^;(Ti — T2) with various combinations of the 
indices suffice), and using all eight permutations of the 
values of the e- variables, we get sixteen different terms at 

order Denoting cos (a% (ri)+A6lfe, (r2)+A6l„„(T3)^ 

by {ij + kl + mn) and so on, and omitting the summation 
and integration signs along with a multiplicative factor 
e:xjp{—3lab/ zoa), they are, 



r f f ^ f If ^ f \ 

— [{ij + kl + mn){Gi^^^i{Ti - T2) - gG(,-fe;(ri - T2) - gG^;„„(T2 - r3) - -G^„^^.(t3 - n) j 

+ (^(U + kl- mn)G{j^^i{Ti - T2) + m.n - u)G{;„„(t2 - T3) - ^(mn -|- ij - kl)G{^„,j{T3 - ti)) 

^{ij -kl + mn)G{^^^i{Ti - T2) + mn - u)(Gf^-fe;(Ti - T2) + G^„^^(r3 - n)) - ^{ij + kl ~ rnn){G{i„^jT2 - T3) 



( 



(B4) 



Keeping track of the signs, using the pairwise symme- 
tries among the indices {{i,j),{k,l),{m,n)), and using 



the property that G/ j.;(ti — T2) is even in space and 
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time, it is straightforward to see that all of the terms 
in Eq. \ti4\ cancel among themselves. Hence, the third 
order cumulant expansion of Eq. |21 does not renormal- 
ize V. Note that, to get this result, we did not have 
to assume any partcular lattice-type; it's generally true 
for all lattices in any dimension. Among other things, 
we have thus confirmed the conjecture made in Ref. [l^ 
about vanishing of the third order term in the corre- 



sponding zero-dimensional problem. Moreover, we have 
established that the same is true in any higher dimension. 
This calculation does not imply, however, that there is 
no order correction to Eq. |S1 As we have seen in 
Eq. 1201 the order term arises already at the second 
order cumulant level, by generating the diagonal Joseph- 
son interaction J in a square lattice. In D = 1, however, 
there is no higher order correction to Eq. [S] 
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